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ABSTRACT 

The exciting detection of a very high degree of linear polarization, P = 80% ± 20%, in the prompt 7-ray emis- 
sion of the recent GRB 021206, provides strong evidence that synchrotron emission is the dominant radiation 
mechanism. Besides this immediate implication, there were also claims that this implies a magnetic field that 
is ordered on large scales within the ejecta, and must therefore be produced at the source, which in turn was 
used as an argument in favor magnetic fields playing an active role in the production of GRB jets. However, an 
alternative explanation was also suggested: a very narrow jet, of opening angle 9j ~ I/7, where 7 > 100 is the 
Lorentz factor during the GRB, viewed slightly outside its edge, at Oj < 8 a \, s < 6j + 1 /j. This explanation also 
works with a magnetic field that is generated in the internal shocks and does not originate at the source. We 
calculate the expected degree of polarization for these two scenarios, and find that it is significantly easier to 
produce P > 50% with an ordered field. More specifically, we obtain P ~ 43 — 61% for an ordered transverse 
magnetic field, B or d, whereas a shock-produced field that is random but fully within the plane of the shock, B±, 
can produce up to P < 38-54% for a single pulse in the GRB light curve, but the integrated emission over 
many pulses (as measured in GRB 021206) is expected to be a factor of ~ 2 lower. A magnetic field normal 
to the shock front, Bu , can produce P ~ 35 - 62% for the emission integrated over many pulses. However, po- 
larization measurements from GRB afterglows suggest a more isotropic configuration for the shock-produced 
field that should reduce P by a factor ~ 2 — 3. Therefore, an ordered magnetic field, B or d, that originates at the 
source, can produce the observed polarization most naturally, while Bi\ is less likely, and B± is the least likely 
of the above. 

Subject headings: gamma rays: bursts — polarization — shock waves — MHD 



1. INTRODUCTION 

The recent detection of a very large linear polarization, 
P = 80% ± 20%, in the prompt 7-ray emission of GRB 
021206 (Coburn & Boggs 2003, hereafter CB), establishes 
synchrotron emission as the dominant radiation mechanism in 
the prompt GRB. As the prompt GRB is believed to arise from 
internal shocks within a relativistic outflow (Rees & Meszaros 
1994; Sari & Piran 1997), it can provide valuable informa- 
tion on the magnetic field structure in the ejecta, and clues 
to the nature of the central source. In a recent paper (Gra- 
not & Konigl 2003, hereafter GK), we suggested that "the 
radiation from the original ejecta, which includes the prompt 
GRB and the emission from the reverse shock (the 'optical 
flash' and 'radio flare') could potentially exhibit a high de- 
gree of polarization (up to ~ 60%) induced by an ordered 
transverse magnetic field advected from the central source". 
This is perfectly consistent with the polarization measured in 
GRB 021206. CB also attributed the polarization in this GRB 
to an ordered magnetic field, and suggested that this implies 
that magnetic fields drive the GRB explosion. A similar inter- 
pretation of this measurement has even been claimed to favor 
Poynting dominated outflows in GRBs (Lyutikov, Periev & 
Blandford 2003). 

However, Waxman (2003) suggested an alternative expla- 
nation: if the GRB outflow is a uniform jet with sharp edges 
and an opening angle 6 ; < 1 /j, then our line of sight is likely 
to be at an angle 8j < 8 b s <8j+l /*y from the jet axis. In this 
case we should see both a bright GRB (as much of the radi- 
ation is still beamed toward us) and a large polarization (e.g. 
Gruzinov 1999; Granot et al. 2002). This scenario does not 
require an ordered field and also works for a magnetic field 
that is generated at the internal shocks (Medvedev & Loeb 
1999). 



There are therefore two feasible explanations for the large 
polarization measured in GRB 021206, where only one of 
them requires a magnetic field ordered on angular scales 
> 1 /j. This undermines the possible theoretical implications 
of an ordered magnetic field in the GRB ejecta. In this Letter 
we critically examine these two scenarios, and estimate their 
ability to explain the high observed polarization. In 50 we 
calculate the polarization from an ordered magnetic field. The 
maximal polarization for a narrow jet with a shock-produced 
magnetic field is calculated in &\ In flwe apply our results 
to GRB 021206 and discuss the conclusions. 

2. AN ORDERED MAGNETIC FIELD 

Here we calculate the linear polarization for synchrotron 
emission from a thin spherical shell with an ordered transverse 
magnetic field, B or( j, moving radially outward with 7> 1. We 
integrate over the emission from the shell at a fixed radius and 
do not follow the different photon arrival times from different 
angles 8 from the line of sight (l.o.s.). This calculation is rele- 
vant to the prompt GRB, the reverse shock (the 'optical flash' 
and 'radio flare') and the afterglow, provided the magnetic 
field is ordered over an angle > I/7 around the l.o.s. 

Following GK, the polarization position angle, measured 
from Z? or d, is given by 8 p = 4>+ arctan(y^-cot0) in the limit 

7^1, where y = (j8) 2 and is the azimuthal angle. We 
have l v = I' v ,{v/v'f, with 1 l' v , oc (^')" Q (sinx') e oc (i/) _Q [l - 
(n 1 ■ B' old ) 2 ]^ 2 , where vjv' w 27/(1 +y), 1 
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1 Here \' is the angle between h' and B', which is also the pitch angle 
between the electron's velocity and B' . For the optically thin part of the 
spectrum that is considered in this work, and as long as the electron energy 
distribution (taking into account electron cooling) is independent of the pitch 
angle \' (which is most natural for a random field, and is also reasonable to 
expect, at least approximately, for an ordered field as well), we find e = 1 + a. 
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Table 1 . Parameter Values for Different PLSs of the Spectrum 
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FIG. 1 . — The polarization P OR | (lower panel) and P ol( ] / Pmax (upper panel) of 
synchrotron emission from an ordered transverse magnetic field, as a function 
of the spectral index a, for e = 1 + a, calculated using Eq. III . 



[(1 — ;y)/(l + y)] 2 cos 2 (j) + sin <j> and n' is the direction in the 
local frame of a photon that reaches the observer. The Stokes 
parameters are given by 
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where a = 3+ a for the instantaneous emission (relevant for the 
afterglow) and a = 2 + a for the time integrated emission (rel- 
evant for the prompt GRB when integrated over a time larger 
than the duration of a single pulse, as in GRB 021206). For 
a uniform jet, the limits of integration should include only re- 
gions within the jet. This is important only if 6* b s + 1/7 > 6j, 
which is expected to be rare for the prompt GRB, but usually 
occurs during the afterglow. When the edge of the jet is at 
y > 3'max ^ a few, the limits of integration may be taken as 

Jo™" d y Jo * d & In this case u = and = -Q/ I =\Q\/ 1 - 
For internal shocks, each pulse in the GRB light curve is from 

a collision between two shells. The emission near the peak of 
the pulse is mainly from 9 < I/7 (y < 1), and may be approxi- 
mated by taking 3'max = 1 . The emission from y > 1 contributes 
mainly to the tail of the pulse. If the latter is included in the 
temporal integration used for measuring P, and is not below 
the background, then we can take y max 3> 1 (the asymptotic 
limit is reached at 3> max > a few). 

Figure [TJ shows /ord/fmax and 2 P or( j as a function of a for 
e = 1 + a (e.g. footnote^, and Table 1 summarizes the results 
for the relevant (optically thin) power law segments (PLSs) of 
the spectrum. 3 

~ Here P m a x = (a+ l)/(a+5/3) = (Peff+l)/(Peff+7/3), where it is useful to 
define p e ff = 2a + 1. For optically thin synchrotron emission, a > — 1 /3, and 
hence P raax > 1/2. This lower limit on P max arises since P = 1/2 is simply the 
low frequency (v <C I'syn) polarization of the synchrotron emission from each 
electron, and therefore P m ax = 1 /2 in PLSs D,E (see Table 1). For PLSs F, G 
and Ft, Pmax is determined by p e ff (=2, p, and p + 1 , respectively), where for 
these PLSs, p c ff is the effective power law index of the electron distribution. 

3 The most relevant case for GRB 021206 is }> max ^> 1 and a = 2 + a, 
for which the approximate formula P or d(a) = 0.016a 4 — 0.052a 3 —0.013a 2 + 
0.335a +0.276, provides a relative accuracy of better than 0.25% for — 1/3 < 



PLS a 



PordlV max — 1) PordtV mas > i) 



D,E 
F 

G 
H 



-1/3 
1/2 
p-l 

2 
P/2 



1/2 
9/13 

g+j 

p+7/3 
p+2 
p+10/3 



0.404 (0.423) 
0.605 (0.623) 

0.605-0.675 
(0.623-0.691) 

0.675-0.726 
(0.691-0.739) 



0.165 (0.306) 
0.435 (0.549) 

0.435-0.563 
(0.549-0.643) 

0.563-0.656 
(0.643-0.709) 



NOTE. — Parameter values for different power law segments (PLSs) of 
the spectrum (that are labeled as in Granot & Sari 2002). Numerical values in 
PLSs G and H are for an electron index 2 < p < 3. The values of P on j without 
(with) parentheses are for a = 2 + a (a = 3 + a) which are appropriate for the 
prompt GRB (afterglow), (see discussion below Eg. [TV 



3. A VERY NARROW JET VIEWED FROM JUST OUTSIDE ITS EDGE 

In this section we calculate the polarization from a narrow 
jet, of opening angle 9j ~ 1 /j, viewed at an angle 9j < f? b s < 
8j+l/-f from its axis. 4 In contrast to ^2j here the magnetic 
field is assumed to be produced at the shock itself, and there- 
fore has symmetry around the direction normal to the shock, 
n s h. Since the more isotropic the magnetic field configuration 
behind the shock, the lower the resulting polarization, we con- 
sider two extreme cases where the field is most anisotropic: 1 . 
a random field that lies strictly within the plane of the shock 
(B = B±, P = P±), 2. a completely ordered field in the direc- 
tion of n sh (B = By, P =P\\). 

Following Ghisellini & Lazzati (1999), we generalize their 
formula so that it would hold for 6> b s > f?,-, 



h JiKTl 0d0I v (ff)P(ff)sin[2il>M] 



P = 



2tt J\8 r 8 0bs \ 



e(Oj-9 ohs )j o a ^ed0U9) + j;^i 



ede 



U0) 



(2) 

where cos^i = (9j-9l bs -6 2 )/29 obs and 9(x) is the Heavi- 
side step function. For Bu we simply have P\\(9) = P max and 
h' ■ B' = n' ■ r = cos 8' w £2 so that l v cx y e/2 /(l +y) 3+a+€ . How- 
ever, for B± we must average over the possible field orienta- 
tions within the plane of the shock: 5 



P±(y) 



f*d</> 



Ay cos 2 <p 



(i+y) 2 



(e-2)/2 



Jo d <f> 



4-vcos 2 <fi 
(l+.v) 2 



'-ft 



(3) 

and l v oc (1 +y) 3 a times the denominator of Eq. 0. For 
e = 2 and 0, P ± (9')/P m . dx = $ = f^gp and -min(y,l/y), 
respectively. Fig. |2]shows —P±(9')/P max for several values 
of e. A larger e implies a larger |P(#')|, as it suppresses l v at 
(n' ■ B 1 ) 2 w 1 where there is a positive contribution to P±(9'). 

a < 3/2. 

4 It is assumed here that the jet has sharp edges, i.e. the emissivity drops 
sharply over an angular interval A6 < 1/7 around 8 = 8j. A smoother edge, 
AO > I/7, would considerably reduce the polarization. 

5 Here P < (P > 0) implies P along (perpendicular to) the plane contain- 
ing w s (, and h'. 
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FIG. 2.— The local polarization, \P±(9')\ = -Pj_(9'), normalized by P max , 
for a magnetic field that is fully tangled within a plane, and emission at an 
angle 6' from the normal to the plane, for e = —0.5, 0, 1,2,3. 



The global polarization from the whole jet is given by 6 
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(4) 



(5) 



(6) 



(7) 



Sobs,/Sj, yj = 



where cost/'i = [(1 -g 2 )y J -y]/(2^ v ^-y), # 

(7^) 2 ,y 1 = (l- ? ) 2 ^,y 2 = (l+^) 2 y > . 

Figures[3]and[4]show P±(q) and P||(g), respectively, for sev- 
eral values of a and yj, using the relation e = 1 + a. For g < 1 , 
\P±\/Pmax 0.2, while rises sharply above q = 1 (the 
larger y^, the sharper the rise), and peaks at q ~ 1 + 1 / y/yj 
(<7 w 1.7—1.8 forjj =1), which is g just above 1 for y^- 1, 
but at q ~ 1 / ^7 ^ 1 f° r X/ ^ 1 ■ The width of the peak 
is ~ 1/ ^fyj, so that the peak is wider (as well as higher) for 
smaller yj. At larger values of q, \P±\ decreases since for 
^obs > (2-3) max(0j , 1 / 7 ) [i.e. q > (2-3) max( 1,1/ yyT)] , the 
jet may be approximated as a point source, and as q increases, 
the emission in the local frame is almost straight backward 
(i.e. h' approaches -n s h and 9' approaches tt), thus suppress- 
ing P±(6') (see Fig. |2j. However, in sharp contrast with 
B±, for Bit even if 9' is only slightly different from w, still 
P\\(9') = Pmax, and Pit approaches P max for q > 2. The tran- 
sition between P»(q > 2) w P max and P||(g = 0) = is very 
gradual for yj <C 1 and very sharp for yj 3> 1 (for which 

—1/2 

the transition occurs at \q — 1| < y,- <C 1). For y^ > 1, 
P||(<7 < 1)/P max < 0.3, which is a little higher than while 
fory 7 - < 1, P {l (q < 1)/P max < 0.6. Fory, > 1 (9j > l/ 7 ) the 
edge of the jet is hardly visible from the interior of the jet and 
P(q < 1 - l/y/yj) « (for both B±_ and B \\ ). 

The above expressions for P±(q,yj, e,a) or P11 (q,yj,e,d) can 
produce afterglow polarization light curves, by using a = 3+ a, 

6 Here, P < (P > 0) means P along (perpendicular to) the direction from 
our l.o. s. to the jet axis. 
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FIG. 3. — The polarization —P±(q)/Pmax for several values of a and yj, 
calculated using Eq. J5J with a = 2 + a and e = 1 +ct. 



e = 1 + a and P m 



ct+l 
a+5/3 



and adding a model for the time 



evolution of 7(f) and 9j(t), which determine q(t) = 9 \ )s /9j(t) 
and yj(t) = [j(t)9j(t)] 2 . One simple model is to assume q(t < 
tj),yj(t > tj) = const, where tj is the jet break time. 7 Note that 
at a fixed observed time, P remains constant within each PLS, 
but changes across spectral breaks. 

4. APPLICATION TO GRB 021206 AND DISCUSSION 

In the prompt GRB the spectral index is usually in the range 

1 /2 < a < 5/4, for which the time integrated polarization (a = 

2 + a, y max 3> 1) from an ordered transverse magnetic field 
(fl ord ) is P ord ~ 43-61% (e.g. Table 1, Fig. □ footnote QJ. 
This is reasonably consistent with the value of P = 80% ± 20% 
that was measured for GRB 021206 (CB). Furthermore, this 
requires a magnetic field that is ordered over angles > I/7 
which can still be -C 9j. 

We now turn to the narrow jet scenario (Waxman 2003). 
For o bs > 9j + 1/7 the observed flux from the GRB drops 
considerably. Therefore, a bright GRB like 021206 requires 
q < 1 + 1/ y/yj. As it is hard to collimate a jet to 9j < l/ 7 , it is 
reasonable to assume 37 > 1 and therefore q<2. Furthermore, 
foryj ^> 1 that is usually inferred from afterglow observations 
(Panaitescu & Kumar 2002), the peak of the polarization is at 
q ~ 1 + 1 / y/y] ~ 1 and has a width Aq ~ 1 / Jyj -C 1 which 
covers a fraction ~ 1 j Jy] <C 1 of the solid angle from which 
the GRB is beamed toward us, and therefore a high polar- 
ization is very unlikely. Hence, we require yj < 2. The fact 
that GRB 021206 was extraordinarily bright, together with the 
correlation found by Frail et al. (2001), might suggest a very 
narrow jet, so that yj < 2 is not so far fetched (Waxman 2003). 
Altogether we expect 1 < yy < 2 and 1 < <? < 2. 

In this parameter range, and for l/2<a<5/4, P± peaks at 
Pi,max ~ (0.55-0.7)P max - 38-54%. However, the Lorentz 

7 Ghisselini & Lazzati (1999) simply assumed 9j,a = const, and implicitly 
assumed e = 2 since they used P(6')JP mw = sin 2 6*7(1 + cos 6 ). However, 
they did not take into account the fact that l' v , oc ((sinx') e ) x ([1 — (P 1 ' 
n') 2 ] e ' 2 ), which effects the polarization light curves. 
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FIG. 4. — The same as Fig.l3*lbut for Pu using Eq. |4}- 



factor of the shocked fluid in the internal shocks is expected 
to vary with A7 ~ 7 between different shell collisions within 
the same GRB. 8 This implies a reasonably large variation in 
y I oc 7 2 , while q = const, so that our l.o.s. will not be near the 
peak of P± for all the pulses in the GRB light curve. Further- 
more, the observed flux at q ~ 1 + 1 / ^/y], where P± peaks, is 
smaller than near the edge of the jet (q w 1), due to relativis- 
tic beaming effects, so that the brightest pulses would tend to 
be relatively weakly polarized, thus further reducing the aver- 
age polarization over the whole GRB. Therefore, while for a 
single pulse in the GRB light curve P± can approach P± imax , 
the average over many pulses (as in GRB 021206) will be 
P±<P± < m*x/2~ 19-27%. 

For By we find P\\ ~ (0.3-0.9)P max - 20-70% for a sin- 
gle pulse, and expect P\\ ~ (0.5-0.8)P max ~ 35-62% for the 
average over many pulses, which is consistent with the value 



measured for GRB 021206. In fact, Bn is an ordered mag- 
netic field, just that unlike Z? or d which was considered in <J3J 
it can in principle be generated at the shock itself, as n s h is 
a preferred direction that is determined locally by the shock 
front. Current models for the production of magnetic fields 
at collisionless relativistic shocks (Medvedev & Loeb 1999) 
suggest B = B± rather than B», However the amplification 
mechanism of the magnetic field and its configuration in rel- 
ativistic shocks is still largely an open question, so that it is 
hard to rule out B sa Bh on purely theoretical grounds. Nev- 
ertheless, it is important to keep in mind that we considered 
two extreme cases for the magnetic field configuration behind 
the shock, in which it is most anisotropic. The relatively low 
values of P < 3% measured in GRB afterglows 9 , compared 
to the expected values of P < 20% (Sari 1999; Ghisellini & 
Lazzati 1999; GK), suggest that the magnetic field created 
behind relativistic shocks is more isotropic than the extreme 
cases we considered, implying P values lower by a factor of 
~ 2-3 (e.g. GK). Therefore, although Pn ~ 35-62%, a more 
isotropic magnetic field configuration that is suggested by af- 
terglow observations would imply 10 P ~ 15 - 30%. 

We therefore conclude that P > 50% is most naturally pro- 
duced by an ordered magnetic field that is carried out with 
the ejecta from the central source (as was recently proposed 
by GK). This is therefore the most likely explanation for the 
value of P = 80% ± 20% (CB) measured in GRB 021206. A 
magnetic field that is generated at the shock itself is less likely 
to produce a sufficiently large polarization. However, if either 
1 ) the systematic uncertainty in the quoted value was for some 
reason underestimated, and P < 20-30% is acceptable, or 2) 
the internal shocks for some reason produce a magnetic field 
much more anisotropic than in the afterglow shock, then Pn 
may still be a viable option. Both points are required in order 
for P± to work well. 
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ences. 



8 If B on j is ordered on angles > l/7min> which are still < 0.01, and B OIL \ 
does not change significantly (i.e. by < 0.5 radians) between the different 
shells, then this should not effect P on j significantly; Pn should also not be 
strongly effected. 

9 except for a possible sharp spike with P S3 10% in the polarization light 



curve of GRB 020405 (Bersier et al. 2003). 

10 The same argument should reduce P± ~ 19 — 27% to P ~ 7 — 13%, 
making it even harder to reconcile with the value measured in GRB 021206. 



REFERENCES 



Bersier, D., et al. 2003, ApJ, 583, L63 

Coburn, W., & Boggs, S.E. 2003, Nature, 423, 415 

Frail, D. A., et al. 2001, ApJ, 562, L55 

Ghisellini, G., & Lazzati, D. 1999, MNRAS, 309, L7 

Granot, J., & Konigl, A. 2003, ApJ, 594, L83 (GK) 

Granot, J., Panaitescu, A., Kumar, P., & Woosley, S. E. 2002, ApJ, 570, L61 
Granot, J., & Sari, R. 2002, ApJ, 568, 820 
Gruzinov, A. 1999, ApJ, 525, L29 

Lyutikov, M., Pariev, V. I., & Blandford, R. D. 2003, ApJ, in press 
Jastro-ph/03054l]5) 



Medvedev, M. V., & Loeb, A. 1999, ApJ, 526, 697 
Panaitescu, A., & Kumar, P. 2002, ApJ, 571, 779 
Rees, M. J., & Meszaros, P. 1994, ApJ, 430, L93 
Sari, R., & Piran, T. 1997, ApJ, 485, 270 
Sari, R. 1999, ApJ, 524, L43 
Waxman, E. 2003, Nature, 423, 388 



